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ABSTRACT 

Strong gravitational lensing is a useful probe of both the intrinsic properties of the lenses and the 
cosmological parameters of the universe. The large number of model parameters and small sample of 
observed lens systems, however, have made it difficult to obtain useful constraints on more than a few 
parameters from lensing statistics. Here we examine how the recent WMAP measurements help improve 
the constraining power of statistics from the radio lensing survey JVAS/CLASS. We find that the absence 
of 61 > 3" lenses in CLASS places an upper bound of /3 < 1.25(1.60) at 68% (95%) CL on the inner 
density profile, p oc r^*^, of cluster-sized halos. Furthermore, the favored power spectrum normalization 
is cTg ^ 0.7 (95% CL). We discuss two possibilities for stronger future constraints: a positive detection 
of at least one large-separation system, and next-generation radio surveys such as LOFAR. 

Subject headings: cosmology: theory - gravitational lensing 



1. INTRODUCTION 

Strong gravitational lensing is a sensitive probe of the 
intrinsic properties of galaxies, clusters, and dark matter 
halos as well as of the underlying cosmological model of 
the universe. For instance, some of the earliest constraints 
on the cosmological constant were obtained from statis- 
tics of gravitational lensing (e.g. Turner 1990; Kochanek 
1996). The uncertainties in the properties of the lenses 
themselves, however, have led some to suggest that lensing 
studies are better used as probes of lens properties rather 
than cosmological parameters (e.g. Cheng & Krauss 2000). 
This view is particularly pertinent in light of the recent 
measurements from the Wilkinson Microwave Anisotropy 
Probe (WMAP) and large-scale structure surveys, which, 
when combined, can now constrain some of the basic cos- 
mological parameters to 5% or better (Spergel et al. 2003 
and references therein). 

Indeed, recent lensing studies have yielded useful con- 
straints on the density profiles and velocity dispersions 
of galaxies (Takahashi & Chiba 2001; Rusin & Ma 2001; 
Oguri, Taruya & Suto 2001; Keeton & Madau 2001; Davis, 
Huterer & Krauss 2002; Chae 2002). Most of these analy- 
ses as well as direct observations (Cohn et al. 2001; Treu & 
Koopmans 2002; Winn, Rusin & Kochanek 2002; Rusin, 
Kochanek & Keeton 2003) indicate that the inner density 
profiles of elliptical galaxies are close to singular isother- 
mal (SIS): p oc r~^ with f3 ~ 2. These analyses have been 
greatly facilitated by the completion of the Cosmic Lens 
All-sky Survey (CLASS; Myers et al. 2003, Browne et al. 
2003), which extended the earlier Jodrell Bank-VLA As- 
trometric Survey (JVAS; King et al. 1999). JVAS/CLASS 
is the biggest survey with a homogeneous sample of lenses, 
with a total of about 16,000 sources and 22 confirmed lens- 
ing events. Among these, a subset of 8958 sources with 13 
observed lenses forms a well-defined subsample suitable for 
statistical analysis (Browne et al. 2003, Table 13). In par- 
ticular, the CLASS statistical subsample finds no lensing 
events with angular separations 9 > 3". An explicit search 
at 6" < 9 < 15" also finds no lensing events (Phillips et 



al. 2002). 

In this paper we study the constraints from strong lens- 
ing statistics in the post- WMAP era. We examine two 
of the currently most uncertain cosmological parameters: 
the equation of state of dark energy w and the power 
spectrum normalization erg, as well as two of the most 
uncertain parameters in current lens models: the inner 
mass density profile of clusters (3 and the mass scale 
that separates galaxy from cluster lenses. A-body sim- 
ulations indicate that the dark matter profile in clusters 
has a roughly universal form with an inner logarithmic 
slope of /3 - 1 to 1.5 (Navarro, Frenk & White 1997; 
Moore et al 1999). However, recent analyses based on 
observations of arcs in clusters (e.g. Sand, Treu & Ellis 
2002) indicate a shallower inner profile. This uncertainty 
has motivated us to treat /3 for clusters as a free param- 
eter in this study. For the galaxy-scale lenses with lower 
masses (M < Mc), we assume the SIS model in accordance 
with the results cited above. This two-population galaxy- 
cluster model is needed to account for the large number 
of observed small-separation lenses and the lack of > 3" 
systems in JVAS/CLASS (Keeton 1998; Porciani & Madau 
2000). 

For the cosmological parameters, we use the results from 
WMAP in conjunction with ACBAR, CBI and 2dF exper- 
iments (Spergel et al. 2003). We use the mean values n = 
0.97, ilMh'^ = 0.134 and f^s/i^ = 0.023 (each measured to 
better than 5%), and marginalize over ^Im — 0.27 ± 0.04 
and as = 0.84 ± 0.04, although we find the strength of 
our final constraints on /3 and Mc to be insensitive to the 
strength of the particular priors used on and erg. We 
consider only flat cosmological models. 

2. LENS PROFILES AND OPTICAL DEPTH 

We write the total optical depth at a given lens mass as 
the sum of the optical depths due to SIS and generalized 
NEW (GNFW) halos: 

t(M) = fsis{M)Ts,s{M) + il-fsisiM))Ta^,,,^{M) (1) 

where fsis is the fraction of halos that are SIS, and 
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1 — fsis is the fraction that are GNFW. In addition to 
these two populations, there is some evidence that objects 
with M ^ 1O"M0 have shallow profiles (Ma 2003, Li & 
Ostriker 2003). Ma (2003), for example, derives a form 
for fsis by requiring the halo mass function and galaxy 
luminosity function to give consistent optical depth pre- 
dictions at subsecond scale. Current lensing surveys are 
limited to 9 ^ 0.3", which cannot constrain small-mass ha- 
los. We therefore use the simplest possible form: fsis = 1 
for M < Mc and otherwise, and leave Mc as a free pa- 
rameter to be constrained by data. 

An SIS lens has a density profile p{r) — a^/{2nGr'^), 
where <Ty is the 1-d velocity dispersion. SIS lenses produce 
an image separation of 26e, where 6e — 'in{ay/c)^Dis/Ds 
is the Einstein radius, and the cross section is cions = 
^{OeDi)'^ = \&^^{(jJc)^{DiDu/Dsf, where Z?„A, and 
Dis are the angular diameter distances to the source, to 
the lens, and between the lens and the source, respectively. 
Using (jy = {■KG^M'^AyirPM/f^y^^, we can then relate the 
cross-section criens to the halo mass M which is needed for 
eq. lHJ. The GNFW profile (Zhao 1996) is given by 
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{r/r,f [l + {r/rs)r 
The mass of a halo is defined to be the virial mass 
M = AirAyirPMf'lir/'^j where r^^^ is the virial radius 
within which the average density is AyirPM and pM is the 
mean matter energy density in the universe. We compute 
Ayir from the spherical-collapse model using the fitting 
formula in Weinberg & Kamionkowski (2003). In prac- 
tice, we replace the scale radius by a concentration pa- 
rameter c{z) = ryir{z)/rs{z), which is well described (for 
/3 = 1) by c(z) = co(l + z)-i(M/M,)-0-i3 with cq = 9, 
where M^, is a typical collapsing mass for that cosmology 
at redshift zero (e.g. Bullock et al. 2001). For (3 ^ 1, we 
compute Co by assuming that the ratio ri/2/ryir is inde- 
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pendent of /?, where 



^1/2 



is the half-mass radius defined 



as M{r < ri/2) = M(r < ryir)/2 (Li and Ostriker 2002). 
The same redshift and mass dependence is used for all p. 

The optical depth for lensing is given by (e.g.. Turner, 
Ostriker & Gott 1984) 



dzi—— (1- 

dzi 



dn 

dM—{M, zi) ai^ns{M, zi)B 
dM 



(3) 

where n(M, zi) is the physical number density of dark ha- 
los at the lens redshift zi, (Tiens(-W, zi) is the lensing cross 
section of a halo of mass M at z;, and B is the magnifi- 
cation bias. We considered accounting for the ellipticity 
of lenses, but found that for the JVAS/CLASS luminosity 
function, ellipticity makes small (^ 10%) correction to the 
optical depth, and therefore we choose to ignore it. This 
is consistent with previous findings that ellipticity mostly 
affects image multiplicities (e.g. Kochanek & Blandford 
1987; WaUington & Narayan 1993). 

To model the number density of lenses, we use the halo 
mass function instead of the galaxy luminosity function 
(for the latter approach to constrain the dark energy equa- 
tion of state, see Chae et al. 2002). The former has the 
advantages that (1) it accounts for all objects - dark and 
luminous - in the universe; (2) it has been accurately cal- 
ibrated by A^-body simulations; and (3) its main depen- 
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Fig. 1. — Number of lenses with angular separation above d from 
the CLASS statistical sample (Browne et al. 2003) and predictions 
from the compound lens model with different values of w), erg, /3 and 
Mc- Halos with M < Mc are assumed to have the SIS profile. 



dence is on now well-measured cosmological parameters 
instead of the more uncertain astrophysical parameters. 
We use the fitting formula of Jenkins et al. (2001) for 
dn/dM but note that their halo mass is independent of 
cosmology and is defined to be the mass enclosed in ra- 
dius Tyir within which the average density is (ISOpcrit); 
instead of [AyirPu) based on the spherical collapse model 
for general cosmology discussed earlier. We use the lat- 
ter to define our halo mass in eq. (3), but for consistency, 
we convert it to Afigo and evaluate the fitting formula for 
dn/dM at Afiso- We find this step to lead to a relatively 
small upward correction to the optical depth. The matter 
power spectrum is an input to dn/dM; we use the for- 
mulae of Eisenstein & Hu (1997) for w = — 1 models and 
those of Ma et al. (1999) for models with w > — 1. We do 
not consider u; < — 1 since the power spectrum for these 
models is uncertain. 

We include a scatter of 0.14 in log^o 2(2) (Bullock et 
al. 2001; note that the published version of this paper 
incorrectly quotes 0.18 for the scatter) for the halo con- 
centration in eq. (3). We find that including this effect is 
very important, as it increases the GNFW optical depth 
by more than an order of magnitude, in rough agreement 
with Chen (2003) and Kuhlen, Keeton & Madau (2003). 

To compute the magnification bias, we note that the 
sources in CLASS are well represented by a power-law lu- 
minosity function, (j){S) — dn/dS oc S~^, with rj ~ 2.1, 
which leads to i? = 4.76 for all SIS lenses (Rusin & 
Tegmark 2001). The magnification bias is more compli- 
cated for GNFW lenses. We use the formula in Oguri et 
al. (2002), which agrees very well with ray-tracing simula- 
tions. 

3. RESULTS 

Before discussing our statistical constraints on (3, Mc, 
as, and w, we first illustrate in Fig. ^ the dependence of 
lensing probability on these parameters by comparing the 
expected number of lenses A'^(> 9) in various models with 
the CLASS result. It shows that the value of Mc and 



3 



10 



10' 





1 1 1 1 1 1 












1 








68% CL 


< 

95% CL 





0.5 



0.75 



1.25 



1.5 



1.75 



Fig. 2.— Shaded regions show the 68% and 95% CL likelihood 
joint constraint on the slope /3 of the density profile of M > Mc halos 
and the mass Mc separating the SIS and GNFW lenses. The results 
have been marginalized over all other parameters using WMAP pri- 
ors. The vertical lines show the upper limits on /3 after further 
marginalization over Me- 



an tipper limit on /3 can be determined very acctirately 
because the CLASS histogram has a sharp cut-off in 9. 
Variations in w, on the other hand, have a fairly small 
effect at ah scales. Models with low as underpredict the 
optical depth a,t 9 ^ 3", while those with both high erg and 
f] overpredict the optical depth for ^ 3". Finally, we find 
that A^(> 9) is insensitive to the value of flM- 

For the statistical tests, we use information from both 
the total lensing optical depth t and the image separation 
distribution dr/dO from JVAS/CLASS. We do not use the 
lens redshift distribution due to selection effects that are 
presumed to be significant in this test. This approach 
is similar to that described in Davis, Huterer & Krauss 
(2002). The combined likelihood for the two tests is 
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where the first term accounts for the total optical depth 
and the second for the angular distribution. Here x = 13 is 
the number of lenses in CLASS, N = 8958t is the number 
of galaxies predicted by the model, and t{w, fjg, (3, M^, Zg) 
is the predicted optical depth computed from eqs. (1) and 
(3). Ti is the optical depth given the source redshift Zg 
for the lens in question. When Zg for a particular lens is 
not available, we set = r, i.e., we set Zg to equal the 
mean redshift of the whole source population, (zg) = 1.27 
(Marlow et al. 2000). 

Fig. 121 shows the joint constraints on /3 and Mc using 
the WMAP priors on Qm and as and the WMAP mean 
values for n, flmh'^ and Q-Bh"^ that are determined very 
accurately. Although we have fixed w = —\ here, we 
have checked that the constraints are insensitive to w by 
repeating the analysis for several values of w and find- 
ing essentially identical results. The absence of 6* ^ 3" 
lenses in CLASS places an upper limit on the optical 




Fig. 3. — 68% and 95% CL likelihood joint constraints on f) and 
(Tg after marginalization over Mc and Ca/- 



depth at these angular scales, which in turn limits how 
steep the inner cluster density profile can be: we find 
(3 < 1.25(1.60) at 68% (95%) CL (vertical lines) after 
marginalization over Mc- The constraints on Mc are also 
strong due to the sudden break in the predicted optical 
depth at the angular separation corresponding to Mc (see 
Fig. 1): we find 1.58 < Mc/ {10^^ H'^Mq) < 2.51 at 68% 
CL and 1.25 < Afc/(10" H-^Mq) < 3.98 at 95% CL after 
marginalization over (3. 

Fig. 121 shows the joint constraints on /? and as using the 
same WMAP results as in Fig.|21and after marginalization 
over Mc- It shows that steep density profiles {(3 ~ 1.5) 
are consistent with the absence of large-separation lenses 
only for a small range of low cg. High as (^ 1), on the 
other hand, is allowed by CLASS data provided the density 
profile is shallow. The lower limit of as > 0.7 at 95% CL 
is insensitive to the value of f3 because lowering as further 
would severely underpredict the cumulative optical depth 
at small image separations. 

We have also attempted to obtain joint constraints on 
as and w using the same WMAP parameters. After 
marginalizing over erg, our likelihood analysis gives w < 
-0.75(-0.55) at 68% (95%) CL. The solid and double- 
dotted curves in Fig. 1, however, show that for two models 
with w = —1 and —0.5 but otherwise identical parame- 
ters, we would expect nearly identical number of lenses at 
9 < 2" and at most a factor of two of difference at larger 9. 
This insensitivity to w is due to the fact that the matter 
power spectrum for different w has identical shapes except 
on very large length scales (see, e.g.. Ma et al. 1999), so 
once the models are normalized to the same cg today, the 
only difference is in the growth rates at higher redshift and 
in for halo masses. (See Sarbu, Rusin & Ma 2001 for 
constraints on w for COBE-normalized power spectrum.) 
The reason we were able to place any constraints on w 
is because the theoretically computed optical depth for a 
fiducial (WMAP-favored) cosmology underpredicts the to- 
tal optical depth, so that various models lie on the tail of 
the likelihood function. Since a more negative w leads to 
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a higher r, the total optical depth part of the likelihood 
function exponentially suppresses models with w > — 1. 
We caution that constraints on w (or any other parame- 
ters to which the test in question is only weakly sensitive) 
are robust only when the theoretically computed optical 
depth roughly agrees with the observation. 

4. DISCUSSION 

We have used the statistics of the recently completed 
JVAS/CLASS radio survey to constrain the parameters in 
the two-population model of halos in the universe. Mo- 
tivated by recent evidence from direct observations, A^- 
body simulations, and semi-analytic arguments, we have 
assumed that all objects with mass less than Mc have the 
SIS profile, while the more massive ones have the GNFW 
profile that scales as in the inner region of the halos. 
The absence of ^ 3" lenses in CLASS enables us to ob- 
tain tight constraints on Mc and an upper limit /3 < 1.6 at 
95% CL. Furthermore, we obtain a constraint of as > 0.7 
(95% CL) on the power spectrum normalization. 

The constraints in our study have come from three ef- 
fects: the total lensing optical depth measured by CLASS, 
the shape of the image separation distribution, and the 
lack oi 9 > 3" lenses. It is therefore interesting to consider 
what happens when at least one large-separation event 
(^ 5") is observed, and when such an event will be ob- 
served given our knowledge of the halo density profiles and 
cosmology. Keeton & Madau (2001) ask a similar question, 
and remark that one such large-separation event should be 
found in surveys such as the Two-Degree Field and Sloan 
Digital Sky Survey (SDSS), or else the cold dark mat- 
ter model would have to be questioned. Fig. ^ however, 
shows that one could easily accommodate an extremely 
small probability of finding a large-separation event if, for 
example, we lived in a universe with erg ~ 0.7. 

This issue may soon be resolved: there is evidence for a 
9 « 13" quadruple lens system in the SDSS (N. Inada et 
al., in preparation). Since this event is not yet a part of 
a controlled survey with significant statistics, it cannot be 
included in our analysis. To test how the constraint on /3 
would change with a positive detection of wide-separation 
lenses, we repeat our analysis by assuming a 14th lens in 
the statistical sample of JVAS/CLASS with 9 = 13" and 
a source at = 2. We find this hypothetical sample to 
favor 1.05 < (3 < 1.35 at 68% CL, while tightening the 
constraint on with its central value unchanged. This 
result illustrates that stronger (lower as well as upper) 
limits on p would be obtained if we had at least one large 
angular separation event in a controlled survey. This is 
simply because models with (3 <, I predict ^ 1 lenses at 
large 9 and would be strongly disfavored. 

These factors motivate us to consider future surveys 
that may have enough statistics to detect large-separation 
lenses. The proposed LOw Frequency ARray (LOFAR; 
www.lofar.org) is expected to find millions of radio sources 
to fluxes below 1/LtJy (Jackson 2002). Although the pro- 
posed angular resolution of ~ 1" is inferior to that of 
CLASS, LOFAR would be a useful survey for detecting 
large-separation systems. Estimates from Jackson (2002) 
indicate that LOFAR would have 17,000 lenses in the 
brighter part of the survey with 30(7 signal-to-noise ra- 
tio, providing a thousandfold increase in lensing statistics 
(if the identity of the lens events can be confirmed). Given 



the small number of lenses in current surveys, even if the 
LOFAR estimates are overly optimistic we can still expect 
that significantly larger statistical lens samples and other 
future wide-field telescopes such as SNAP (snap.lbl.gov) 
and LSST (www.dmtelescope.org/dark_home.html) will 
greatly improve our understanding of the populations of 
objects in the universe. 
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